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ABSTRACT 

We have calculated the temperature and density structure of the hot postshock plasma 
in magnetically confined accretion flows, including the gravitational potential. This 
avoids the inconsistency of previous calculations which assume that the height of the 
shock is negligible. We assume a stratified accretion column with 1-d flow along the 
symmetry axis. We find that the calculations predict a lower shock temperature than 
previous calculations, with a flatter temperature profile with height. We have revised 
previous determinations of the masses of the white dwarf primary stars and find that 
for higher mass white dwarfs there is a general reduction in derived masses when the 
gravitational potential is included. This is because the spectrum from such flows is 
harder than that of previous prescriptions at intermediate energies. 

Key words: accretion, accretion disks - methods: data analysis - novae, cataclysmic 
variables - stars: fundamental parameters - white dwarfs - X-rays: stars 



1 INTRODUCTION 



Magnetic Cataclysmic Variables (mCVs) are strong emitters 
of hard X-rays. These are produced in the postshock region 
immediately above the accreting white dwarf. In this post- 
shock region there is a strong temperature gradient between 
the shock front (> 10 8 K) and the surface of the white dwarf. 
By correctly modelling the properties of the emission region 
it is possible to determine fundamental proporties of the 
mCV, such as the mass of the white dwarf. For an overview 
of the modelling the hard X-ray emission see Cropper et al. 
1999, while Cropper (1990) and Warner (1995) provide more 
general reviews of mCVs. 

Wu (1994) (also Wu, Chanmugam & Shaviv 1994) de- 
scribed a closed integral solution for the temperature and 
density profile of this postshock flow. The treatment in- 
cluded cyclotron cooling as well as thermal bremsstrahlung: 
cyclotron cooling is important for those mCVs with strongly 
magnetic (> 10 MG) white dwarfs. The advantage of this 
scheme is that the temperature and density profiles can be 
calculated analytically (with the exception of a final sum- 
mation). Thus multi-temperature spectra can be calculated 
sufficiently rapidly to permit spectral fits to X-ray data. 
Previously, an analytical solution existed only for the pure 
bremsstrahlung case (Aizu 1973). The assumptions for the 
boundary conditions in Wu (1994) are the same as those 
in Aizu (1973): these are a cold preshock flow (at free-fall 
preshock velocities) followed by a strong shock and then a 
hot postshock flow cooling onto a cold white dwarf surface 



(reaching zero temperature and velocity at the base of the 
flow). As in Aizu (1973), the height of the shock was as- 
sumed to be negligible so that the effects due to gravity 
were neglected. 

Cropper, Ramsay & Wu (1998) and Ramsay et al. 
(1998) used the Wu (1994) formulation to determine the 
masses of the white dwarfs in those mCVs observed with 
Ginga. The masses they derived are at the higher end of 
expectations, and in the case of XY Ari, are significantly 
higher than the best determinations obtained by other 
means (Ramsay et al. 1998). This indicates that there is 
more hard X-ray flux than expected, or less soft X-ray flux 
as a result of significantly more complex forms of absorption 
than assumed in Cropper et al. (1998). The latter has been 
explored for BY Cam in Done & Magdziarz (1998). The ef- 
fect of absorption in mCVs is a topic of its own, and we do 
not pursue it here. Issues of the emission from a structured 
accretion region can be dealt with to first order by the sum- 
mation of a sufficiently large number of local models with 
different accretion rates and magnetic fields. This is because 
the emission from the postshock flow is optically thin in the 
continuum down to almost the white dwarf surface. There- 
fore, here we take the path of exploring the effect of reducing 
the restrictiveness of the assumptions in the emission model 
itself. 

Several options present themselves at this stage. These 
include improvement of the lower boundary condition to 
match the atmosphere of the white dwarf more appropri- 
ately, separate treatment of the proton and electron popu- 
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lations and more detailed treatment of the physics within 
the shock itself. It should be noted that there have been a 
number of calculations for the postshock radial flows onto 
white dwarfs, for example Imamura & Durisen (1983) (and 
references therein) and Woelk & Beuermann (1996). Some 
of these include 2-fluid effects, Compton cooling and the ef- 
fect of a gravitational potential. These are all important. 
The aim of those calculations has been to determine the 
postshock temperature and density, and predict spectra and 
X-ray light curves, but they are unsuitable for the iterative 
model fitting of X-ray data. This differs from the approach 
adopted by Wu (1994) and Cropper et al. (1998) which at- 
tempts to extract information from the X-ray spectra and 
which requires a formulation that can be computed suffi- 
ciently rapidly for that purpose. Here we continue along the 
path of improvements to that technique by addressing the 
elimination of a negligible shock height assumption. In so 
doing we explore and elucidate clearly for the first time the 
effects of including a radially varying gravitational acceler- 
ation on the postshock flow structure. 

The justification for assuming that the shock height 
has negligible effect is given, for example, in Frank, King 
& Raine (1992). This is adequate for the first generation 
of mass estimates as in Cropper et al. (1998): however, in 
those cases where the cyclotron cooling is insignificant (the 
IPs) and where the mass of the white dwarf is larger than 
~ O.8M , application of the Aizu (1973) or Wu (1994) for- 
mulae results in very significant shock heights for typical 
specific accretion rates (0.3 7? WD for a 1.0 Mq white dwarf 
accreting at 1.0 g s _1 cm -2 ). This is inconsistent with the 
assumptions used by Aizu (1973), Wu (1994) and Cropper et 
al. (1998) (where the gravitational acceleration is ignored) 
and with Woelk & Beuermann (1996) (who use a constant 
gravitational acceleration) . Since the temperature in the im- 
mediate postshock region is the pre-shock velocity at that 
height (divided by a factor 4 because of the strong shock 
jump condition) this implies that the shock temperature will 
be significantly lower than that calculated from the free fall 
velocity at the surface as used when assuming a negligible 
shock height. For this paper, we have therefore augmented 
the treatment in Wu (1994) to include the effect of the vari- 
ation of the gravitational acceleration within the postshock 
flow. 



2 FORMULATION 

We derive the following set of 1-d steady state conserva- 
tion equations taking into account the gravitational poten- 
tial (see Appendix for details): 



d 2 -GMp 



dP „dv 
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With the ideal gas law, P/p = kT/p,mH, this set of equa- 
tions is closed. Here x is the spatial coordinate, p is the 
density, v is the flow velocity, P is the pressure, 7 is the adi- 
abatic index, T is the temperature, /1 is the mean molecular 



mass and rriH is the mass of a hydrogen atom. A is the cool- 
ing term, which includes both bremsstrahlung and cyclotron 
radiation (Wu 1994), G is the constant of gravitation and 
M is the mass of the white dwarf. We have included in Ap- 
pendix A the derivation of the form of the A term (used but 
not shown in Wu 1994) for the reader to assess the treatment 
of the cyclotron cooling. Note that we limit ourselves strictly 
to a 1-d flow: we do not for this work consider spherically 
symmetric accretion or accretion in a dipolar field geometry, 
as the form of (1) is particularly useful for our purposes. 

The temperature and density structure of the postshock 
flow can be computed from the coupled pair of nonlinear first 
order ODEs derived in Appendix A. Rapid increases in the 
density at base of the flow make it more advantageous to 
integrate using the velocity v as the independent variable: 
(A9) should therefore be inverted and then the dependence 
for the composite variable £ (see Appendix A) can be written 
in terms of the product of (A8) with the inverse of (A9). We 
therefore have 
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Here — C is the mass transfer rate, v a — 4v s is the free- 
fall velocity at the height of the shock, A is the coeffi- 
cient for bremsstrahlung radiation, a — 3.85 and [3 = 2 
are the coefficients corresponding to cyclotron radiation in 
Wu, Chanmugam & Shaviv (1994) and e s is the ratio of the 
bremsstrahlung to cyclotron cooling time at the shock, as in 
Wu (1994). 

Such a coupled system requires two boundary values. 
Unfortunately this is not an initial value problem: the initial 
value x — R WD at v — is known, but the variable £ is not - 
it is constrained at the shock front by the requirement that 
the pressure immediately behind the shock P a = 3Cv a /4. 
Moreover, this is a floating boundary condition as the freefall 
velocity v a at the height of the shock is itself a result of the 
computation. 

We have also combined the coupled pair (4 and 5) into a 
single second order nonlinear ODE in x(v) (A14). Formally 
this now becomes an initial value problem as at x — -R WD , 
v — (as above) and ^ = 00. However because of poles 
in (A 14) neither of these are accessible. Stepsizes must be 
chosen appropriately in order to maintain the stability of 
the integration and this route therefore appears to provide 
no computational advantage at this stage. 



3 IMPLEMENTATION 

We have used a Runge-Kutta-Merson method in a shooting 
and matching technique implemented in the NAG D02HBF 
routine (Numerical Algorithms Group 1995) for floating 
boundary conditions to solve the coupled pair of equations 

4 and 5. Initial guesses for the parameters of the routine are 
derived from the zero-gravity case. Convergence is generally 
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Figure 1. Upper plot: The temperature profile for the case in- 
cluding the gravitational potential (solid line) compared to the 
standard zero gravity profile (dotted line) when cyclotron cooling 
is included (30 MG) and negligible (0 MG). The vertical lines are 
at the position of the shock front. Lower plot: as above but for 
the electron number density profiles. See text for details. 



Figure 2. Upper panel: The photon spectrum from the Aizu 
profile (above), and, displaced downwards by a factor of 10 for 
clarity, that for the case when the effects of gravity are included. 
The parameters in Figure 1 are assumed here, and the system 
distance is 100 pc. Lower plot: the ratio of the Aizu spectrum to 
the spectrum from the varying gravitational potential. 



reached within a few iterations, thus permitting a rapid cal- 
culation of the postshock structure. The additional compu- 
tation time overhead by comparison with the closed integral 
form in Wu (1994) is negligible. 

We show in Figure 1 the temperature and density profile 
of the postshock flow for the zero gravity case and our re- 
vised case including the gravitational potential. We assume 
a 1M© white dwarf accreting at 1 g s _1 cm 2 , correspond- 
ing to an accretion rate of ~ 4 x 10 15 g s _1 over a fraction 
/ = 0.001 of the white dwarf, typical for mCVs. The MG, 
zero gravity case is the standard Aizu profile, while the 30 
MG zero gravity case is the same as the Wu (1994) profile. 

The effect of treating the gravitational potential more 
appropriately is significant, especially in the case of the pure 
bremsstrahlung cooling (0 MG field). The shock temper- 
ature is lower than that calculated from the Aizu profile, 
because the pre-shock velocities at this height are less than 
those at the surface of the white dwarf. However, there is 
significant heating of the postshock flow itself through the 
release of gravitational potential energy, so the temperature 
in the subsequent flow is higher than that calculated from 



the Aizu profile. This leads to a flatter profile. In the case 
where we include cyclotron cooling from a 30 MG field, the 
correction is not so significant, as the height of the shock is 
in any case lower. 

In Figure 2 we show the spectra for the pure 
bremsstrahlung case. The ratio plot (lower panel) indicates 
that the new profile produces a harder spectrum with fewer 
lines, at least for the 0.1 to 10 keV range shown here. This 
means that fits to X-ray data using the new profile will, 
somewhat counterintuitively, produce masses for the white 
dwarf which are lower than those determined using the Aizu 
(1973) or Wu (1994) profiles (Cropper et al. 1998). 

At higher energies, the Aizu spectrum becomes harder, 
as the postshock flow immediately behind the shock is at a 
higher temperature. However, because for IMq white dwarfs 
the shock temperature is ~ 50 keV, this is somewhat higher 
than the typical range sampled by typical datasets available 
for these systems at present. 
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4 IMPROVED FITS TO DATA 

Cropper et al. (1998) and Ramsay et al. (1998) derived 
masses for the accreting magnetic white dwarfs in the Ginga 
archives. Some of these exceeded 1.0 Mq . This indicates 
that the shock height in these particular systems is likely to 
be significant, especially for the IPs in the sample. There- 
fore their white dwarf mass determinations can be improved 
significantly by taking into account the revised temperature 
and density profiles calculated above. 

We have proceeded as in Cropper et al. (1998) and Ram- 
say et al. (1998), so that a direct comparison can be made 
with the earlier mass derivations. The results are shown in 
Table 1 for the former paper and in Table 2 for the lat- 
ter. An intermediate level of refinement is also shown in 
the tables: the earlier papers assumed a pure Hydrogen 
plasma (p, = 0.5) for calculating the temperature and den- 
sity profiles of (but not the spectra from) the postshock 
flow. Tables 1 and 2 therefore show the original mass deter- 
minations, the revised determinations with a cosmic plasma 
where fi — 0.615, and then the present calculations including 
the gravitational potential and /i = 0.615. 

In addition to cold interstellar absorption, Cropper et 
al. (1998) and Ramsay et al. (1998) used a partially ionised 
absorber to account for the absorption within the system of 
the preshock flow. In the presence of density inhomogeneities 
in the accretion flow, and arc-shaped accretion footprints, 
the treatment of the absorption is a complex matter (Done 
& Magdziarz 1998, Rainger 1998, in preparation). The final 
ingredients for the absorption model would include ionised 
absorption for the finely divided flow and a range of partial 
covering fractions for the density inhomogeneities, all inte- 
grated over a range of pathlengths dependent on the viewing 
angle to the flow. As noted earlier, this is beyond the scope 
of this work. However we have included in Table 1 a final 
column which shows the fits including the gravitational po- 
tential and n = 0.615 but now assuming a partial covering 
of cold material, rather than the partially ionised absorber. 
This may be useful for comparison with other workers, as the 
partially ionised absorber model is not generally available. 

As expected from Figures 1 and 2, the results in Ta- 
bles 1 and 2 indicate reduced masses for the massive white 
dwarfs in these systems when the gravitational potential is 
included in the shock structure. In the case of Table 1 the av- 
erage (unweighted) reduction in mass of those systems above 
1.0 Mq in Cropper et al. (1998) is approximately 0.1 Mq . 
For those below 1.0 Mq the masses are slightly (~ 0.03Mq ) 
increased. There is little difference between the masses as- 
suming jj, = 0.5 and n — 0.615 in the zero gravity case. There 
is some scatter in the trends because the fitting routine op- 
timises parameters such as mass transfer rate in addition 
to the white dwarf mass, and in some cases new global op- 
tima are found. The masses derived using the partial fraction 
covering model for the absorption are slightly lower still: for 
those systems previously above 1.0 Mq the average reduc- 
tion in mass is a further ~ 0.12 Mq , while that for the less 
massive systems is unchanged. 

In the case of Table 2, we also see reductions ~ 0.1 
Mq as the effect of the gravitational potential is included in 
the mass determinations for XY Ari. 



5 CONCLUSIONS 

The results in Tables 1 and 2 show a reduction in the mean 
of the derived masses for the white dwarf when the effects of 
gravity are included. The reductions are greatest for more 
massive white dwarfs, low specific accretion rates and low 
magnetic fields, because the shock height is greater under 
these conditions. The determinations for lower mass white 
dwarfs remain largely unaffected. 

Comparing the revised determinations with those from 
other methods, we note that agreement in the case of XY 
Ari is significantly improved: for the RXTE data the best 
mass and 2<r errors are in close agreement with those deter- 
mined from the eclipse duration in Ramsay et al. (1998). The 
Ginga and Asca determinations remain significantly higher, 
however: Ramsay et al. (1998) explored some of the possible 
explanations for this. In the case of the systems in Table 
1 for which other determinations exist (see also Cropper et 
al. 1998) we note close agreement in the case of the lower 
mass systems EX Hya as determined from the X-ray line 
data (Fujimoto & Ishida 1997, Ishida 1999, Ezuka & Ishida 
1999), and with AO Psc and TX Col (Hellier et al. 1996, 
Ezuka & Ishida 1999). In the case of FO Aqr, where the 
mass is higher at ~ 1 Mq , there is also close agreement with 
that determined from Asca data (Ezuka & Ishida 1999). The 
Ginga data for AE Aqr and QQ Vul are relatively poorer, 
and do not constrain the masses strongly: both are never- 
theless consistent with those determined elsewhere (Welsh, 
Home & Comer 1994, Mukai & Charles 1987). Only in the 
case of AM Her is the mass we derive still significantly higher 
than those from other techniques (Gansicke et al. 1998, and 
to a lesser extent, Mukai & Charles 1987). However the ef- 
fect of the absorption model is particularly strong, so that if 
a partial covering model is assumed, even here the mass of 
0.85 M Q is not controversially high (Mukai & Charles 1987). 

We note that any technique for determining the mass of 
the white dwarf will have both systematic biases built into 
the method and uncertainties arising from poor constraints 
in some of the input parameters, such as inclination for the 
kinematically derived masses, or perhaps instrumental ef- 
fects as in the case of XY Ari (Ramsay et al. 1998). 

Our results indicate the importance of the absorption 
model on the mass determinations. The effect can be min- 
imised by excluding the softest photons from the fits to the 
data, especially in the case of CCD data with sufficient spec- 
tral resolution. The resulting models can be used retrospec- 
tively to isolate the absorption at lower energies. Clearly, 
with shock temperatures in excess of 40 keV for 1 Mq white 
dwarfs, X-ray spectra with sufficient data quality extending 
to approximately equivalent energies are beneficial for the 
extraction of mass information. 
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System 


,u=0.5+IA 


^=0.615+IA 


/i=0.615+G+IA 


/i=0.615+G+PC 




M Q(x 2 red ) 


M @(x 2 red ) 


M @(x 2 red ) 




M e(x 2 ed ) 


AM Her 


1.22 (1.43: 1.15-1.30) 


1.12 (1.28: 1.06-1.20) 


1.06 (1.35: 


0.97-1.16) 


0.85 (0.61: 0.76-0.90) 


EF Eri 


0.88 (1.00: 0.82-0.94) 


0.81 (1.07: 0.76-1.15) 


0.80 (1.02: 


0.75-0.85) 


0.80 (0.91: 0.68-1.17) 


BY Cam 


1.08 (1.53: 1.00-1.25) 


1.27 (1.43: >0.63) 


1.18 (1.38: 


0.76-1.27) 


0.98 (1.16: >0.65) 


V834 Ccn 


0.5 (0.88: 0.25-1.00) 


0.5 (0.82: <1.3) 


0.64 (0.89: 


<1.15) 


0.54 (0.95: <1.25) 


QQ Vul 


1.13 (1.40: 0.95-1.30) 


1.30 (1.32: >0.7) 


0.95 (1.27: 


>0.6) 


0.60 (1.18: 0.4-1.4) 


EX Hya 


0.52 (0.65: 0.42-0.60) 


0.45 (0.59: 0.32-0.52) 


0.50 (0.59: 


0.44-0.56) 


0.46 (0.56: 0.42-0.50) 


AO Psc 


0.40 (0.93: <0.80) 


0.45 (0.90: <0.62) 


0.36 (0.89: 


<0.65) 


0.56 (0.94: 0.36-0.72) 


FO Aqr 


1.22 (0.83: >1.00) 


1.12 (0.77: >0.89) 


1.07 (0.73: 


0.91-1.22) 


0.92 (0.73: 0.58-1.22) 


TV Col 


1.20 (0.84: >0.95) 


1.22 (0.93: 1.03-1.36) 


1.21 (0.82: 


>1.06) 


1.3 (1.21: >0.9) 


BG CMi 


1.25 (0.62: 1.05-1.39) 


1.36 (0.93: >1.18) 


1.19 (0.95: 


>0.95) 


1.09 (0.58: >0.94) 


TX Col 


0.55 (1.34: 0.44-0.64) 


0.48 (1.28: 0.41-0.54) 


0.48 (1.23: 


0.37-0.53) 


0.48 (1.26: 0.39-0.54) 


PQ Gem 


1.35 (1.08: > 1.05) 


1.32 (1.04: >1.15) 


1.21 (1.06: 


>1.08) 


1.29 (1.04: >1.17) 


AE Aqr 


0.3 (0.76: <1.15) 


0.62 (0.68: 0.4-1.15) 


0.6 (0.69: 


<1.0) 


0.6 (0.66: <1.0) 



Table 1. The best fit to the mass of the white dwarf in the magnetic CVs which were observed using Ginga using variants of the 
stratified accretion column model of Cropper ct al. (1998). The upper panel shows the results of the synchronous (or very near) systems 
(the Polars), while the lower panel shows the results of the asynchronous systems (the Intermediate Polars). The column with the mean 
molecular mass of the plasma, fi=0.5, corresponds to that assumed in Cropper et al. (1998). (Some of these new values differ marginally 
to those quoted there since in this paper the accretion shock has been divided into 100 vertical elements as opposed to 50 there). The 
remaining columns assume fi=0.615 with the last 2 columns including the gravitational term described in this paper (+G). IA refers to 
an ionised absorber while PC refers to partial covering of cold absorber. The range in mass is the 90 per cent confidence interval. In all 
cases a cold absorber (for interstellar absorption) and a fluoresence line fixed at 6.4keV have been included. 



Satellite 


fi = 0.5+IA 


,u=0.615+iA 


/i=0.615+IA+G 




M @ (x 2 red ) 


M © (x 2 red ) 


M e (x 2 red ) 


RXTE 








(1) 


1.03 (1.04: 0.94-1.11) 


0.98 (1.04: 0.85-1.03) 


0.89 (1.04: 0.84-1.00) 


(2) 


0.89 (0.76: 0.83-1.02) 


0.83 (0.76: 0.75-0.92) 


0.82 (0.76: 0.74-0.91) 


ASCA 








(GIS2) 


1.05 (0.88: 0.81-1.32) 


1.01 (0.88: 0.75-1.25) 


0.96 (0.89 


0.72-1.16) 


(GIS3) 


1.38 (1 


06: >1.18) 


1.37 (1.07: >1.12) 


1.25 (1 


08 


1.05-1.36) 


(GIS2&3) 


1.25 (0 


97: 1.08-1.38) 


1.23 (0.97: 0.98-1.36) 


1.15 (0 


98 


0.96-1.25) 


(SIS0) 


1.26 (1 


05: >0.96) 


1.26 (1.05: >0.86) 


0.98 (1 


07 


0.78-1.32) 


(SIS1) 


1.40 (1 


11: >1.27) 


1.40 (1.10: >1.23) 


1.31 (1 


10 


>1.07) 


(SIS0&1) 


1.40 (1 


05: >1.24) 


1.40 (1.05: >1.17) 


1.19 (1 


08 


1.00-1.36) 


(GIS&SIS) 


1.33 (1 


03: >1.22) 


1.31 (1.03: >1.13) 


1.15 (1 


05 


1.02-1.27) 


Ginga 










1.28 (0.46: >1.17) 


1.23 (0.46: 1.08-1.38) 


1.19 (1.19: 1.06-1.32) 


Hcllicr (1) 




0.91-1.29 




Hellier (2) 




0.74-1.14 




Refining 


1 & 2 


0.78-1.03 





Table 2. A comparision of mass estimates using XY Ari data using the models described in Table 1. No Iron fluorescence component 
at 6.4keV has been added. The range in mass is for the 90 per cent confidence level. The data in the first column have previously been 
reported in Ramsay et al. (1998). 



REFERENCES 

Aizu, K., 1973. Prog. Thcor. Phys., 49, 1184 

Chanmugam, G., Wu, K., Courtney, M.W., Barrett, P. E., 1989, 
ApJS, 71, 323 

Cropper, M., 1990, Sp. Sci. Rev., 54, 195 

Cropper, M., Ramsay, C, Wu, K., 1998, MNRAS, 293, 222 

Cropper, M., Wu, K., Ramsay, G., 1999, in Hellier, C, Mukai, K., 
eds. ASP Conf. Ser., Proc. 2nd Workshop on Magnetic Cat- 
aclysmic Variables, Annapolis, July 1998, Astron. Soc. Pac, 
San Francisco, in press 

Done, C, Magdziarz, P., 1998, MNRAS, 298, 737 



Ezuka, H., Ishida, M., 1999, preprint 

Frank, J., King, A. R., Raine, D., 1992, Accretion Power in As- 
trophysics, Cambridge Univ. Press, Cambridge, 2nd edition 

Fujimoto, R., Ishida, M., 1996, ApJ, 474, 774 

Gansickc, B., Hoard, D. W., Bcucrmann, K., Sion, E. M., Szkody, 
P., 1998, A&A, 338, 933 

Hellier, C, Mukai, K., Ishida, M., Fujimoto, R., 1996, MNRAS, 
280, 877 

Imamura, J. N., Durisen, R. H., 1983, ApJ, 268, 291 

Ishida, M., 1999, in Hellier, C, Mukai, K., eds. ASP Conf. Ser., 
Proc. 2nd Workshop on Magnetic Cataclysmic Variables, An- 



© 0000 RAS, MNRAS 000, 000-000 



6 M. Cropper et al. 



napolis, July 1998, Astron. Soc. Pac, San Francisco, in press 
Mukai, K., Charles, P., 1987, MNRAS, 226, 209 
Ramsay, G., Cropper, M. S., Hellier, C, Wu, K., 1998, MNRAS, 

297, 1269 
Wada, T., Shimizu A., Suzuki M., Kato M., Hoshi R., 1980, Prog. 

Theo. Phys., 64, 1986 
Warner, B., 1995. Cataclysmic Variable Stars, Cambridge Univ. 

Press, Cambridge 
Welsh, W. F., Home, K., Comer, R., 1994, MNRAS, 257, 649 
Woclk, U., Beuermann, K., 1996, A&A, 306, 232 
Wu, K., 1994. Proc. Astron. Soc. Australia, 11, 61 
Wu, K., Chanmugam, C, Shaviv, C, 1994, ApJ, 426, 664 



APPENDIX A: HYDRODYNAMIC EQUATIONS 

The time dependent mass continuity, momentum and energy 
equations are (see e.g. Frank, King & Raine 1992): 

dp 



dt 



+ V • (pv) = o 



P^+pv-Vv = -VP + f 
at 



^Pv 2 + pe) + V • [(^pv 2 + P e + P)v] = 
f • v - V • F rad . 
Here we ignore heat conduction and fluid viscosity. 



(Al) 
(A2) 

(A3) 



Al The Radiative Cooling 

By Gauss' Theorem, we can replace the volume integral of 
the divergence of the radiation flux V • F ra d in (A3) by a 
surface integral 



d A x V ■ F r . 



ds -F r 



Suppose the radiative flux F rat j consists of of two com- 
ponents: an optically thin bremsstrahlung term Fbr and an 
optically thick cyclotron term F cyc , thus, 

F ra d = Fbr + F c yc • 

Since bremsstrahlung emission is optically thin and 
isotropic, the surface integral of its flux at any arbitrary 
closed surface S" is 47r times the volume integral of the 
bremsstrahlung emissivity in a volume V' enclosed by the 
surface. Therefore, we have 



■ F cyc /47r to the 
cyclotron emissivity j cyc . As the volume integral of the di- 
vergence of the flux is the same as the integration of the 
flux over the surface enclosing the volume, we need only to 
evaluate the flux at the surface of the cylinder. Cyclotron 
emission is strongly beamed such that most power prop- 
agates preferentially in the direction perpendicular to the 
magnetic field, the vector products of ds • F cyc at the top 
surface 5*1 and the bottom surface S3 of the cylinder are 
small in comparison with that at the side surface 5*2 (this 
is not equivalent to assuming that the fluxes at the surfaces 
Si and 53 are small). As an approximation we neglect their 
contribution to the surface integral. Then we have 



d 3 x V- 



F c 



ds -F c 



where V c is the total volume of the cylinder. 

Suppose the spectrum of the cyclotron emission is opti- 
cally thick up to a frequency v*{— cj„/2tt) after which it is 
optically thin. For parameters typical of the accretion shocks 
in mCVs the optically thin cyclotron intensity falls rapidly 
with frequency as ~ z/~ 8 (see Chanmugam et al. 1989). As 
a first approximation we can neglect the contribution of the 
optically thin cyclotron emission to the total cyclotron cool- 
ing process determining the structure of the postshock ac- 
cretion flow. Thus 



ds-F c 



*2 



dx I duB R - I (v) 



where 73 RT = 2kTu 2 /c 2 is the Rayleigh- Jeans intensity, and 
x s is the shock height and D is diameter of the cylinder 
(which is the accretion column) respectively. 

If we divide the accretion column in layers of height dx, 
the cyclotron luminosity, which is the total energy loss due 
to cyclotron emission, is the sum of the contribution of these 
layers: 



dx^^ 
dx 



Thus, we obtain an effective local cyclotron cooling term 



_ riLcyc _ n kTuoj 
cyc ~ dx _7r 127T 2 C 2 ' 

Now from Wada et al. (1980) 



d A x (47rj br ) = / d x A 



ds ■ F br = 



where Ab r is the energy loss per unit volume due to 
bremsstrahlung emission, which is the bremsstrahlung cool- 
ing function. Since the volume V' is arbitrarily chosen, we 
have in all space 



r e - 

l_10 4 _ 


0.05 


r T 1 
_io 8 kJ 



w.(x) ~9.87w c 



where the cyclotron frequency u) c = eB/m c c and the di- 
mensionless plasma parameter = 2nn c (x)D/B (normally 
written A). Here B is the magnetic field and n c is the elec- 
tron mumber density. For constant B, 



V-F br 



A b 



V 



where A is the constant for bremsstrahlung emission, P is 
the pressure and p is the density. 

To evaluate the volume integral of the divergence of the 
cyclotron flux, we consider the postshock emission region 
as a cylinder with its symmetry axis parallel to the mag- 
netic field. Because cyclotron emission is not optically thin 



We define two quantities 
3 1 

fcyc = -(lie + Th)kT— 

•^ Jlcyc 

and 

3 1 

£ br = -(n e + n p )kT- — 

l Abr 
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for the cyclotron and bremsstrahlung cooling respectively 
(rii is the ion number density) . In terms of these two quan- 
tities, the total effective cooling term is 

A = A br + A,,, = A,,, ( I + -^- 

with 

£br 



We assign a proportionality e(x) — tbi/t cyc and scale it in 
terms of the pressures and densities at the shock, e s . With 
the ideal gas law, we eliminate the temperature and obtain 

3.85 

e(x) 



ibi 



tc 



(0 



Thus, we arrive at the effective cooling term determining the 
dynamics of the flow as that used in Wu (1994) 



A = Ab 



l + e s TT ^ 



(A4) 



A2 Steady State 1-D Hydrodynamic Equations 
including Gravity 

For a stationary state, the time-derivative in the hy- 
drodyanamic equations are zero. As we consider a 1- 
dimensional flow channelled along a cyclinder parallel to the 
magnetic field, (Al) to (A3) can be reduced to 

d 



dx {pv) = ° 

dv dP _ 
dx dx 



(A5) 
(A6) 
(A7) 



The symbols are as defined in Section 2, and / is the force 
term, e is the internal energy of the gas. As the cross sec- 
tion of the accretion column is small in comparison with the 
white dwarf radius, it is adequate to let 



/ = pg 



pGM 



and in (A5), (A8) and (A9) we recover the hydrodynamic 
equations of Wu (1994), but with the additional gravita- 
tional force term in the momentum equation. 



A3 Coupled First Order Form 

We now substitute the variable £ = v + (P/ pv) in place 
of P in the hydrodynamic equations. Thus the momentum 
equation (A8) becomes 



d , „, -GMp 
Using (1) we obtain 



d£ 

pv — 
dx 


-GMp 

X 2 


and thus 


d£ _ 


-GM 



dx x 2 v 

Substituting for £ in the energy equation (A7) and not- 
ing from (Al) that pv =constant= C (where —C is the mass 
transfer rate M), 

„d£ ^dv dv , 

Thus 



dv 

dx 



O x* 

7(1 ~v)-v 



(All) 



Now substituting for Abr in (A4) we have 

0' 



A = Ap 2 



thus 



l + e s 



A& + P 



(£-u)V 



:i a 



„ a +P 



(A12) 



where we have used zero pre-shock pressure (so that across 
the shock from (2) p^v 2 = p B v 2 + P B and p a w a = p B v B ) to 
eliminate P s and p s and symbols are defined in Section 2. 
Combining (A10) to (A12) we obtain equations (4) and (5) 
in Section 2. 



dv 

dx 



d .(p V 2 + P)^-, l - 



Noting from (A5) that 
d 2 d dv 

we obtain for (A6) 

GM 
dx^ ' ' x 2 

Substituting for / from (A6) in (A7) yields 

d , 1 2n A 2 7-,n dv 

V dx- {P£ -2 PV) + { 2 PV +P£ + P) Tx=- K - 

Using e — — —37, and again (A5), we obtain 

v dP 7P dv 

+ - L -r— =-A 



7 — 1 dx 7 — 1 dx 
and thus 
dP „ cfo , , , , 

v te + ^ P Tx = -^- 1)K 



(A8) 



(A9) 



A4 Second Order Form 

We now proceed to the 2nd order form in v. We eliminate £ 
by substituting further w — (£/v) — 1 and 

A«+0 1 



., a +P ' 



Equation (4) becomes 



dx 
dv 



v("/w - 



-(7- WvMl + <^ 5 



+ 



(A13) 



For e' a = (no cyclotron) this is a quadratic in yfw and two 
roots can be found analytically by rearranging (A13) and ap- 
plying the standard formula for quadratic roots. Otherwise 
this is a quintic with no general analytic form for the roots, 
which then have to be found numerically for the particular 
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Differentiating (A13) again we have 



<Px_ 
dv 2 



-(7- 



^ACJw~ GM 
1 Z— + ^r- 



( , yw - 



-(7-1) 



AC 



dw 

1 ~d~v~~ 

dw 



2jw~ 



2^/w dv 
GM /2dx l\"|-i 
x 2 v \x dv v J J J 



,.ACJw GM 
-(7-1) 5 H— 5- 



where, from (5) 

dw _ GM dx w + 1 
dv x 2 v 2 dv v 



(A14) 



(A15) 



Thus by substituting for y/w from the roots of (A13), and 
^ from (A15) into (A14) we obtain a 2nd order nonlinear 
ODE in x and v alone. 
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